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Abstract 

Mathematical models of biological growth commonly attempt to distinguish deformation due to growth 
from that due to mechanical stresses through a hypothesised multiplicative decomposition of the deformation 
gradient. Here we demonstrate that this hypothesis is fundamentally incompatible with shear- resistance and 
thus cannot accurately describe growing solids. Shifting the focus away from the kinematics of growth to the 
mechanical energy of the growing object enables us to propose an "energy-deformation decomposition" which 
accurately captures the influence of growth on mechanical energy. We provide a proof and computational 
verification of this for tissues with crystalline structure. Our arguments also apply to tissues with network 
structure. Due to the general nature of these results they apply to a wide range of models for growing 
systems. 
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1. Introduction 

1.1. Background 

Biological growth — of cells, tissues, organs and organisms — leads to morphological change as well as 
mechanical stresses such as tension in arteries (see, for example, Holzapfcl and Ogdcn (2010)) and plant 
stems (for example, Goriely and Neukirch (2006), Vandiver and Goriely (2008)). These play an important 
role in biological function. Modeling biological growth and the accompanying mechanical stresses is of 
increasing interest in the biological, continuum mechanics and mathematical communities as it becomes 
abundantly clear that mechanical stresses are not only passive responses to growth, but also feed back 
to influence morphological development as well as biochemical pathways; see, for example, Lintilhac and 
Vesecky (1984), Lynch and Lintilhac (1997), Huang (2004). 

Continuum models for growth typically draw from approaches first developed in the context of non- 
biological continua, for example, plasticity theory or mixture theory; see Goriely et al. (2006), Ambrosi 
et al. (2011) for reviews. A very popular current approach is based on a multiplicative decomposition of the 
deformation gradient. The idea, introduced by Rodriguez et al. (1994), is as follows: Suppose that a map if 
describes the deformation of a body from a reference configuration to a current configuration. In the absence 
of growth, the elastic energy density of the current configuration is a function W of the deformation gradient 
F :— Dip. If, however, growth has also contributed to the deformation, the mechanical energy is a function 
of only that part of the deformation resulting in elastic stresses. Analogously to a standard approach in 
plasticity theory (Lee (1969)), the approach of multiplicative decomposition posits that the gradient is a 
product 

F = AG (la) 
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of tensors G, arising from growth, and A, arising from elastic deformation. (The decomposition (la) is 
understood to hold point-wise in space-time.) There is not, in general, a deformation of the body with 
gradient G; that is, growth may not be compatible with any actual deformation. The elastic deformation 
associated with A restores compatibility (that is, the gradient nature) of F. The elastic energy is a function 
only of A: 

W = W{FG'^). (lb) 

1.2. Main results 

From (lb), the multiplicative decomposition can be understood as saying that energetically growth is 
a motion (in its domain, by G) of the energy density W. Here we present a (in our view) compelling 
argument that this cannot be true in the presence of shear resistance (see Definition 1.1 below) and non- 
dilational growth (defined in Section 1.4 below). Moreover we also propose a way to overcome this: that in 
addition to the deformation gradient being decomposed (multiplicatively) , the energy density be decomposed 
(additively). 

This is done in Theorems 4.1 and 4.12, and Corollary 5.1. These focus on growing systems with a 
crystalline structure; we explain the reason for this below. However our insights and arguments are general 
and we discuss growing networks and growing continua in Section 6. We highlight three features of our 
approach: 

First, we study continua via discrete systems: We begin with discrete systems because they provide 
a context in which the concept of growth can be clearly and unambiguously formulated, and where the 
interaction of growth with deformation can be rigorously derived. To understand growth on a continuum 
level, we compare the continuum limits of the initial and the grown system. (However, in some instances, 
discrete systems themselves are appropriate models for growth with the discrete elements corresponding, for 
example, to biological cells, see Odell et al. (1981), Weliky and Oster (1990), Muhoz ct al. (2010).) 

Second, we shift attention away from the kinematics of growth, cf., (la), to the energetics of growth, 
cf., (lb). We submit that not only is this a better conceptual approach to growth but also facilitates the 
use of more sophisticated mathematical tools such as discrete-to-continuum limits and, more generally, the 
tools of the calculus of variations. In addition it enables a closer integration of the mechanics of growth with 
the biochemistry of growth, and thus permits the development of more holistic models of biological growth. 

Third, we recognise that there are different kinds of growth processes, two of which we discuss in Section 2 
below. 

1.3. Organisation of the paper 

We begin, in Section 2, by considering growahle springs. These are one-dimensional elastic objects which 
are able to change their rest-length by a non-elastic process, namely growth. There is more than one way 
in which this can occur; we present two ways, which we name replication (Section 2.1) and recombination 
(Section 2.2). Of these replication, which involves mass transfer, more naturally corresponds to an intuitive 
understanding of growth. However, recombination, which is a constant-mass process, may be thought of 
as a one-dimensional conceptualisation of the biological process sometimes called "remodelling", see, for 
example, Taber (1995, 2001), Ambrosi et al. (2011), Menzel and Kuhl (2012). Recombination is also simpler 
to analyse mathematically so we focus on it in this paper. However, as a reader who follows our arguments 
can easily see, our claim that the multiplicative decomposition cannot adequately describe growing systems 
is true for generic growth processes including replication. Similarly, while we develop our alternate approach 
of "energy-deformation decomposition" in the context of recombination, the insight underlying it applies to 
generic growth processes (albeit the resulting formulation might be less elegant for replication). 

In Section 3.1 we introduce (node-spring) lattices, which provide our model of a biological system. We 
do this, not because lattices are good models for biological systems but rather because they present the 
simplest context in which we can communicate our insights and develop our arguments. It would have 
been biologically more natural to use (node-spring) networks but we judged that the resulting need for 
more detailed mathematical analysis would have obscured rather than clarified the essential features. After 
this, in Section 3.2, we introduce a one-dimensional example which serves as a concrete context in which to 
introduce the questions that would occupy our attention. 
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Next we come to the heart of the paper: In Section 4 we present our main results (in the context of 
homogeneous lattices growing homogeneously) and arguments (which are much broader in scope). Section 5 
extends these results to inhomogeneous situations by demonstrating that the energy-deformation decompo- 
sition is stable under homogenisation, both of growth and of elasticity. While Sections 4 and 5 considered 
continuum limits of lattices, in Section 6 we briefly touch upon the application of our insights to networks 
and continua. 

Notation and definitions 

Z is the set of integers, N the set of positive integers and is the set of non-negative real numbers. We 

set 

:=Z^n [0,7V]-°, 
dZ^ ■=Z° nd[0,Nf, 

where D is the dimension of the space and N E N. We denote the standard Euclidean basis in by 
{ei, 62, . . . en}- For brevity we set 

Giihj ■ — i 7 
^iij±k • — 6^ lb Cj lb 6/i; 

for i,j, k ~ {1, 2, 3}. (•, •) is the standard Euclidean inner product on M^. For 6* G M, 

\^sm^ cost' J ^ ' 

We denote the identity matrix (in any dimension) by /. An isotropic matrix is one that is a multiple of the 
identity. A dilation is a deformation whose gradient is isotropic. 
The cardinality of a set S* is 1 5*1 . 

We use the term shear, and the related phrases vanishes on a shear and shear-resisting, in a specific 
sense in this paper: 

Definition 1.1 (Shears). (A non-trivial deformation gradient) F e M^^^ \ SO{D), is a shear if 

\\-\\oF=\\-\\ 

on a basis ofM.^. That is, there exists a basis {vi, ■ ■ ■ ,u_d} o/M^ such that 

llFv.ll = yi^l,...,D. 

Now let W: M^^-° ^ M with W{I) ^MW . 

1. W vanishes on a shear if there exists a shear F such that W{F) = W{I). 

2. W is shear-resisting (or resists shears) if there is no shear F for which W{F) = W{I). 

Note that if e M^^-° is a shear then so is 50(2) F S0{2). 
2. Growable springs 

Our discrete systems are built of growable springs which are one-dimensional objects which are charac- 
terised by a rest-length, £, and a current length, e. In addition to elastic deformation, which changes e but 
not £, the springs can also grow, i.e., change £. 

The elastic energy W of such a spring depends on both £ and e. The nature of this dependence itself 
depends on the nature of the growth process being considered. In a moment we consider two kinds of ^owth 
processes, replication and recombination. However elasticity motivates the following constraints on W: 
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Figure 1: Replication 



We require that when e — i, W attains its minimum, which for convenience we set to be 0. Further we 
ask that this state be a stable elastic equilibrium: 

w{e,e) = o, (2a) 

^Jv{£,e)^0, (2b) 

^wii,e)^o. (2c) 

More generally we require that W{£, •) be convex. 

Note that W is the elastic energy stored in the spring. Any energy expended in growth can be accounted 
for by the inclusion of a growth energy which can be added to W to yield the total energy. 

Next we consider two kinds of growth processes which we name replication and recombination. Although 
an analysis of mass-transfer is outside the scope of this work it is natural to think of recombination as a 
mass-preserving process and replication as changing the mass of the system. However here we are concerned, 
not with the underlying mechanisms which accomplish growth, but only with the interaction of growth (of 
either kind) with elasticity. 

2.1. Replication 

We introduce replication through an example; see Figure 1. 

Consider two systems. The first, the initial system consists of a single spring of rest length £. Its elastic 
energy is Wi{-) = W{£, •). The other, the grown system^ consists of two such springs attached in series. Its 
(total) rest-length is 2£. Let its elastic energy be Wg = W{2£, •). 

Let the grown system be elastically deformed to a (total) length of 2e. Since W{£, ■) is convex, an 
energetically optimal state of the grown system is when each spring deforms to length e. That is, 

Wg{2e) = W{2£, 2e) = 2W{£, e) = 2Wi{e). 

In other words, the grown system is energetically equivalent to two copies of the initial system. We call the 
corresponding growth process (i.e., the process that transforms the initial system into the grown system) a 
replication. More generally: 

Definition 2.1 (Replication). A replication is a growth process for which the elastic energy is a (positively) 
1-homogeneous function: For all •q,e,£ > 0, 

W{ri£,r^e) =riw{£,e). (3a) 
Let W{-) := W{1, ■). From the 1-homogenity of W , 

W{£,e)^£W{j). (3b) 

Note that the mass of the grown system is £ times the mass of the initial system. 
Remark 2.2. Note that 1-homogenity requires that the minimum elastic energy be 0: 

minW{£,e) = W{£,£) = £W{1). 

£,e 
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Since this is true for all £ > we conclude that W{1) = 0. 

Moreover convexity ofW{£,-) implies convexity ofW. From (2), 

^Wii,i) ^0=^ W"il) > 0. 

Note that 

^^W{£,£) = W{1)~W'{1)^0, 

Thus e = £ is a stable equilibrium for the growth process. 
Example 2.3 (Hookcan springs). A Hookean spring has energy 

where k{£) is the spring constant of a spring of rest-length £. Considering (as above) such springs joined in 
series leads to the well-known scaling k{£) = £~^k{l). Thus 

W{£,e)=£\k{l){j-l)\ 



which corresponds to 



W{-)^\k{l){--l)\ 



Remark 2.4 (Linearisation). Let £o be the initial length of a spring. For small growth £ = £o -\- 7^o (with 
7 <C 1 j and small deformation e = £o + e^o ( with e ^ 1) we have 

Wi£o + l£o, io + e£o) = {£o + iQW « £oW{l + e - 7) « ^^W" {l){e - ^f. 

2.2. Recombination 

As with replication, we introduce recombination through an example; see Figure 2. 

Consider two systems. The first, the initial system consists of two springs in parallel, each of rest length 
£. Its elastic energy is W^{-) ^ 2W {£,■). 

The other, the grown system also consists of two such springs, but attached in series. Its (total) rest- 
length is 2£. Let its elastic energy be Wg{-) = W{2£, •). 

Let the grown system be elastically deformed to a (total) length of 2e. As before, the optimal configu- 
ration is for each spring to deform to length e. Thus, 

Wg{2e) = W{2£, 2e) = 2W{£, e) = W,{e). 

In other words, the grown system is energetically equivalent to the initial system. We call the corresponding 
growth process a recombination. More generally: 

Definition 2.5 (Recombination). A recombination is a growth process for which the elastic energy is a 
(positively) 0-homogeneous function: For all ri,e,£ > 0, 

W{Tj£,Tje) = Wi£,e). (4a) 

Let W{-) :— W{1, •). From the 0-homogenity of W , 

W{£,e) = W{j). (4b) 



%JWA-m^VWv-9-WA-%-VWv- 

i 



%JWA-m^VWv-9-WA-%-VWv- 

i 



%JWA-m^vwv-9-WA-m-vwv- 
0-m-#Jww#^vwv#^/ww 

i 

; 

Figure 2: Recombinations 



Remark 2.6. As before convexity ofW{£,-) implies convexity ofW, and, from (2), 

W{e,£) = =^ 14^(1) = 0, 

^w{e,i) = o =^ w'{i) = o, 



de 

W{£,£) ^ =^ W"{1) ^ 



Note that 

Thus e = £ is a stable equilibrium for the growth process. 

Example 2.7 (Hookean springs). Consider a rearrageable spring system which consists of N Hookean 
springs, each with spring constant k{£o) and rest-length £„. When the rest-length of the system is £ < N£o, 
the system consists of parallel systems of j- springs in series. When the total length of the system is e 
the energy of the system is 

■ N p 

W{£,e)^-k{£,)£l{-^lf 

which corresponds to 

W{-) = ^k{£o){--l)\ 



Remark 2.8 (Linearisation). Let £o be the initial length of a spring. For small growth £ — £o + 7^0 (with 
7 <C 1 j and small deformation e = £0 + eio ( with e <^1) we have 

W{1, + 74, 4 + e4) = W (^q^) « W^(l + e - 7) « \w"[r){e - 7)^. 

2.3. Generic growth 

Replication and recombination are special cases of growth processes that can be represented by a (posi- 
tively) p-homogeneous energy: 

W{t,e)^£mCj) 

with recombination corresponding to p = and replication to p = 1. 

When the rest-length of such a spring changes from £i to £g the elastic energy density changes from 

W.{-)^ilW{j) 

to 

w,{-)^tiw{-). 

Thus 

where g = > is the growth of the spring. Note that we allow g G (0, 1). 

3. Preliminaries 

3.1. Lattices 

While it excludes, for example, triangular and hexagonal unit cells, the following restricted definition of 
lattices suffices for our purposes: 

Definition 3.1 (Lattices). A lattice is a triple {lP,C,C) where 

1. the connectivity C C is finite and satisfies 

V eC ^ -V iC (5) 

and 

2. the rest-length L is a function 

X — y ^ C or ?/ — a; e C} (0, 00). 

We think of a lattice as a node-spring system where the nodes are points in TP and x,y G are 
connected by a spring iJfx — yGC ory — xGC, in which case, for brevity, we call {x,y) € Z^ x Z^ a spring. 
The rest length of spring {x,y) is given by C{x,y). Where C is irrelevant we speak of the lattice (Z^,C). 

For V £ C we say that v is nearest-neighbour if \\v\\ ~ 1 and next-nearest neighbour if \\v\\ = \/2- 

(The reason for the restriction (5) will become clear in Definition 3.3 below.) The simplest lattices are 
those that are homogeneous: 

Definition 3.2 (Homogeneous lattices). A lattice is homogeneous if C is translation-invariant, that is, 
C{x,y) is a function only of x — y. In this case we view C as a function on CU —C and write L(x — y) for 
C{x,y). 

A critical property of a lattice is its order: 

7 



Definition 3.3 (Order of a lattice). The order of a lattice (Z^,C) is the smallest integer K eN for which 
there exists a partition Ci, . . . ,Ck of C such that each Ci is a set of linearly-independent vectors in M^. 
(Recall that a partition of a set S is a collection of subsets Si, ... , Sk of S which are mutually disjoint and 
whose union is S.) Note that K '-jj . 

The central object of our study is the continuum energy density of a lattice: 

Definition 3.4 (Continuum energy density of a lattice). Let {'E^,C,C) be a lattice. Then its continuum 
energy density W(c,£) • K"^^"^ — >■ M is given by 

W(c.)(f^):= lim 4 inf Vw^fM^^^^M), (g) 

"''N x — y£C 

if it exists. 

Since they concern elasticity and not growth per se, we do not here entertain questions regarding the 
existence and form of W(c,£) but refer the reader to E and Ming (2007). 

Henceforth we assume that the lattices we consider possess continuum energy densities. In addition we 
assume that the homogeneous lattices we consider satisfy the Cauchy-Born rule, see E and Ming (2007). It 
follows that the continuum energy density of a homogeneous lattice is of the form 



W(cx)(F) = Ew^(S)- (7) 

vec ^ ' 



We shall refer to such a continuum energy density as a Cauchy-Born continuum energy density. 

We also prefer to avoid questions that concern elastic homogenisation and not growth per se. To do this 
it is convenient to consider homogeneous representatives of inhomogeneous lattices: 

Definition 3.5 (Homogeneous representative of an inhomogeneous lattice). Let (Z^,C, C) be an (inhomoge- 
neous) lattice and {Z^,C,C) a homogeneous lattice. {Z^,C,C) is a homogeneous representative of {T^^ ,C,C) 
if their continuum energies agree, that is, if 

(These energies are given by (6) and (7) above.) 

The following Lemma can be understood as saying that the continuum energy density of a homogeneous 
lattice describes a solid only if its order is greater than one: 

Lemma 3.6. The Cauchy-Born continuum energy density of a homogeneous lattice is shear-resisting only 
if the order of the lattice is greater than one. 

Proof. Consider a homogeneous lattice {Z^ ,C) of order 1 with C = {vi, V2, ■ . . , v\c\}. Then its Cauchy-Born 
continuum energy density, W, is of the form 

I'^l /II F 



z— 1 ^ 

for some Li, i = 1, . . . ,\C\. We show that W vanishes on a shear. 

Since the lattice is of order 1, |C| ^ D. If |C| < I? we extend C to a basis C := {vi,V2, . . . ,vd} of . 
Pick Rc G SO{D) such that {v\,RcV2) 7^ {vi,V2). Let Fc G K^^^ be the linear operator that satisfies 

Fcvi = vi 

FcVi^RcVi, i = 2,...,D. 

Note that Fc ^ SO{D) since {Fcvi,FcV2) = (ui,i?c«2) 7^ ("i, W2). On the other hand, clearly || -HoF = || -H 
on C. Thus Fc is a shear and W{Fc) = W(/). It follows that W is not shear resisting. □ 
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Figure 3: The one-dimensional lattice (Z,{ei}). 



3.2. An introductory example: The one- dimensional lattice (Z,{ei}) 

We begin with a one-dimensional discrete system that consists of -I- 1 nodes located on Zjv with 
neighbouring nodes connected by springs of rest-length L, see Figure 3. This is the initial (discrete) system. 
Now let the springs grow as follows: The rest-length of the spring connecting nodes located at j — 1 and 
j becomes {g{-^) — L where g: [0, 1] — ^ M+ is a specified increasing function. This is the grown 

(discrete) system. 

When the displacement of the boundary of the system is specified, Uq = and uat = FN, F > 0, the 
energy of the systems is obtained by a minimisation over the displacements of the interior nodes: 



N 



Wl^(F):^ min V , 

uo=0 J = i 

UN=FN 



N 



W^''{F):= min V g(^)-5(^) W. ,,-i,w) 

UN=FN 

where W is the energy density of the springs and the subscripts i and g denote the initial and grown systems. 

It is an easy exercise to show that (with some coercivity conditions on W) these discrete energies converge 
to a continuum limit: 

W^(F)^ lim -Wl'iF)^ min W ( dx 

N~^oo N u: [0,1]^R Jo \ L J 

«(0)=0,m(1)=_F 



Wg{F)= lim ^WN{F)= min [\p{x)W (^^^\ dx 



where G = g' . 

Thus the grown continuum elastic energy-density of the system Wg is related to the initial continuum 
elastic energy-density Wi through: 

Wgi-) = GPW,{-G-') (8) 

where Wi = W. Thus when p = (that is, for recombination) this system satisfies the multiplicative- 
decomposition first introduced in the context of morphoelasticity in Rodriguez et al. (1994). 

This illustrates, in the context of a simple example, the central question of this paper: How is the 
continuum elastic energy-density of the grown system related to that of the initial system? In particular is 
it related through quantities, such as G in this example, that can be thought of as continuum descriptors of 
the growth of the system? 

4. Homogeneous lattices recombining homogeneously 

4.1. The two-dimensional lattice (Z^, {ei, 62, ei±2}) 

Next we explore the simplest two-dimensional lattice that supports elasticity, (Z^, {ei, 62, 61^2}) i the two- 
dimensional square lattice with nodes located at and with nearest-neighbour and next-nearest-neighbour 
interactions, see Figure 4. 
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Figure 4: The two-dimensional lattice Z2 



Let the initial (i.e., at some instant during growth) rest-lcngths of the springs be Li, L2, ii+2 and Li_2 
for the horizontal, vertical, right-diagonal and left-diagonal springs respectively; note that these are arbitrary 
but fixed. At some later instant let the growth (due to recombination) in the springs be gi, 32, 51+2 and 
gi-2 for the horizontal, vertical, right-diagonal and left-diagonal springs respectively. The following result 
relates the continuum limits of the initial and grown systems: 

Theorem 4.1 (Homogenenously-recombined (Z^, {ei, 62, ei-i-2}))- Let the initial and homogeneously-recombined 
(Z^, {ei, 62, 61^2}) lattices have Cauchy-Born continuum energy- densities Wi and Wg respectively. Then: 

1. There exist Wi, W2 '■ M^^^ — > M (independent of gi, g2, gi+2 o-nd gi-2) and Gi, G2 G M^^^ (independent 
of Li, L2, Li+2 and Li^2) such that VF e M^^^, 

W^{F)^Wl{F) + W2{F), (9a) 
Wg{F) = Wi{FG^^) + W2{FG2^). (9b) 

2. // there exists G G M^''^ such that VF e M2x2^ 

Wg{F) = W,{FG-') (10) 

then, in fact the lattice has dilated under recombination and G can be picked to be a multiple of the 
identity. 

Clearly Gi and G2 describe the growth of the system at the continuum level and Item (1) introduces an 
energy- deformation decomposition that relates the initial and grown systems. Item (2) on the other hand 
shows that this is generically optimal and that a multipicative-decomposition of the deformation gradient 
can be expected to hold only when the growth is isotropic, see Example 4.7 below. In the sequel when 
we refer to "energy-deformation decomposition" we mean a result of the form (9) and when we refer to 
"multiplicative decomposition" we mean a result of the form (10). 

Before we present the formal proof we first explain the central insight (Remark 4.2), present explicit 
expressions for Wi, W2, Gi, G2 (Lemma 4.3) and correct what appears to be a popular misunderstanding 
in Remark 4.5. 

Remark 4.2 (Insight behind the proof). Since the springs lie along Ca, a G {1,2, 1 ± 2} (in the reference 
configuration), the elastic energy density of the lattice depends on a deformation F (only) through the 
quantities ll-FeaH, ^ {IjS, 1±2}. Had these vectors formed a basis for the growth of these springs could 
have been represented by a linear operator on M? . But this generically not the case and thus, in general, 
710 G G M^^^ can geometrically represent the growth of the system. However from the four vectors we can 
form two sets of basis and thus the growth of the springs can be represented by two linear operators on M^, 
namely, Gi,G2 G M^^^. Item (1) then follows from an elementary calculation. 

Suppose, however, that the growth ga, a G {1,2, 1 ± 2} of the springs can in fact be represented by a 
linear operator G G M^^^, that is, suppose that there exists aG E M^^^ such that g^ — Ge^, a G {1, 2, 1±2}. 
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Another elementary calculation shows that, nevertheless, (10) does not hold except when G \\ I. The reason 
for this is that the energy density of a lattice with springs along Gca, a € {1,2, 1 ± 2} is not related to the 
energy density of a lattice with springs along Ca, a G {1, 2, 1 ± 2} through a linear change of variables except 
when G \\ I . 

Thus analogies to crystal plasticity ( where the deformation gradient is decomposed into a plastic-part and 
an elastic-part) are misleading since crystal plasticity does not change the (global) metric on the lattice hut 
growth does. 

The next Lemma shows that Wi and W2 can be picked to be purely elastic (i.e., independent of growth); 
and conversely Gi and G2 can be picked to be purely geometric (i.e., independent of elasticity): 

Lemma 4.3 (Non-uniqueness of decomposition). One choice of {Wi,Gi) and (W2,G2) is 

W.iF)^w(m)^w(\\^\. G.^r^; (lla) 



(lib) 



Li J \ L2 J ' ' \0 92 

V Li+2 J V ^1-2 / ' V 31-2/ " 

^ 1 / 51+2 + .91-2 gi+2 - gi-2 
2 V51+2 - .91-2 91+2 + 91-2 

where W is the energy density of the springs. 

Two other genuinely different (i.e., not obtainable by addition of constants) choices exist: 

W,iF) = w(^^)+w(^^'^), G,^(3^ 5i+2-.9r\ ^^^^^ 



and 



W2{F)^[^)+w(^^^), G2^(j\' (12b) 



Li J \ Li+2 / ' \^ .91+2 

2V 

L2 J ' " \ Li_2 7' " \92-91-2 92 



Li J V ^1-2 /' VO 51-2 

W2{F) = W[^^\+w(^^^^\, G2^( °'). (13b) 



Li+2 ) \9i+2 - 92 92 

Remark 4.4. Each of the six energies introduced in Lemma 4-3, in (11), (12) and (13) vanishes on a 
one-parameter family of shears: For (11), Wi and W2 vanish on 



30(2) il ^°'^),eeM and SO{2) { 



1— COS 1+cos 6 

f2 V2 
;in 6 sin 9 

■J2 V2 



respectively; for (12), Wi and W2 vanish on 



respectively; and for (13), Wi and W2 vanish on 
respectively. 
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Remark 4.5 (Residual-stress and shear-resistance). Note that Theorem 4.I holds regardless of whether the 
system, either initially or when grown, is residually stressed or not. In particular (10) holds for dialatory 
growth even in the presence of residual stress and does not hold for non-dialatory growth even in the absence 
of residual- stress. This shows that, contrary to what appears to he popular opinion, the essential issue is 
not the build-up of residual stress under growth. Rather it is the interplay between shear-resistance and 
growth. Indeed the multiplicative-decomposition is valid for a homogeneously-growing fluid as Lemma 4-13 
below illustrates. (See also Example J^.l and §-^.i.-/ below.) 

We are now ready to prove Theorem 4.1 and Lemma 4.3: 

Proof. The Cauchy-Born continuum energy densities are 

W.iF)^ wO^), (14a) 



Qe{l,2,l+2,l-2} 



W,{F)^ Y: Wr^g-'). (14b) 

Qe{l,2,l+2,l-2} 



Now let Wi, W2, Gi and G2 be as in (11). From (14a) it is immediate that (9a) is satisfied. Thus, to prove 
Item (1) it suffices to verify (9b): For F e 
Since G^^Ca = 5q ^ea, a € {1, 2}, 



p2x2. 



W{'^gA+w(^^g^^). (15a) 



Similarly, since G2 — Qa^^on a e {1 ± 2}, 

W2iFG^^) = W (W'^BS^) + ^ 

\ Fi+2 J V -^1-2 

= W (^.^ri) + W {}^9-A) ■ (15b) 

Equations (14b) and (15) immediately yield (9b). 

Similarly, when Wi, W2, Gi and G2 are as in (12), (9a) is immediate and to verify (9b) it suffices to 
verify that 

Gi^Sa^ g^^Sa, a e {1,1 + 2}, 
G2^ea= g^^Ca,, a e {2, 1-2}. 



This in turn immediately follows from 



-1 -1 -1 
9i 9i+2 - 9i 

3r+2 

r^-^ _ ( .9i-2 
^2 ^ — 1 —1 — 1 

\92 -9l-2 92 



Finally, when Wi, W2, Gi and G2 are as in (13), (9a) is again immediate and to verify (9b) its suffices 
to verify that 

Gi^ea^g'^Ca, a e {1,1 -2}, 
G^^e^ = <?-ie„, a e {2, 1 + 2}. 
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This in turn immediately follows from 

^2 ^ — 1 —1 —1 

V.91+2 - .92 92 

We now turn to Item (2). Suppose there exists G G M^""^ such that (10) holds VF e R'^''^. That is, 

Qe{l, 24+2, 1-2} \ " / ae{l,2,l+2,l-2} ^ " / 

Recall that G is independent of Li, L2, ii+2 and Li_2- Since each i„, a G {1,2, 1±2} is arbitrary (though 
fixed) and appears only in one term it follows that the equality holds term- by-term: VF e M^^^, 

\\FG-'eJ = \\Feo,\\g-\ a G {1,2,1 ±2} (16) 

With the choice F = I, setting H = FG^^, this implies that 

^ ■\H{ei+e2)f = \\Hei\\^ + 2{Hei,He2) + \\He2\\^ = \+2{Hei,He2) ^ 



Thus, 



31+2 31 .92 

^ = ||i/(ei - 62)11' = WHe^W^ - 2{HeuHe2) + HffesH' = ^ - 2{Hei,He2) + \. 
9i-2 9i 92 

111 1 



gf 9l gl+2 9I-2 



(17a) 



With the choice F = G, (16) implies that 

25?+2 = l|G(ei + e2)|p - ||Gei||2 + 2(Gei, Ge2) + \\Ge2f - gf + 2(Gei, Ge2) + gl 
2gl_2 = ||G(ei - e^)^ = ||Gei||2 - 2(Gei,Ge2) + ||Ge2|P ^ gl - 2(Gei, Ge2) + 52- 

Thus, 

91+91^91+2+9^2- (17b) 
With the choice F = ^"J, setting H = FG-\ (16) implies, 

^1±A ^ \\H{ei+e2W = \\Hei\\^ + 2{Hei,He2) + \\He2\\^ = 2 + 2{Hei,He2), 

91+2 

+ ^ _ g^^|j2 ^ ||^g^||2 _ 2(iJei, i/e2) + ||i?e2f = 2 - 2(iJei, iJes). 



5L2 



Thus, 



(3?+.9i) -^ + -^ =4. (17c) 

\i/l+2 i/l-2/ 

An easy calculation shows that (17) implies that gi ~ g2 = gi+2 = 31-2 =: g (say). In other words, the 
lattice has dilated by g and, from (16) we can pick G = gl. □ 

Remark 4.6. The proof shows why decomposition into two energies suffices: As Remark 4-4 shows each 
of Wi and W2 in (9) vanishes on a shear. This sufficiently decouples the energy that each "growth mode" 
does not influence the energy of the other "deformation mode", i.e., the deformation that the other energy 
function can see. Thus Gi does not influence W2 and vice versa. 

On the other hand that T4^i(-G^^) and VK2('G^^) are evaluated at the same F couples these energies 
together. Indeed the theorem proves that this provides all the coupling that is needed. 
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Remarks 4.2, 4.5 and 4.6 can be illustrated by a one-dimensional example: 



Example 4.7. Here we present a (geometrically) one- dimensional example where the multiplicative-decomposition 
fails hut an energy- deformation decomposition continues to hold. Consider a one- dimensional discrete sys- 
tem where both the nearest neighbours and the next-nearest neighbours are connected by springs. This is the 
initial system. 

Let the nearest-neighbour .springs grow by gi and the next-nearest neighbour-springs grow by 52 • This is 
the grown system. 

Then, as is easy to see, the continuum energy densities are related by 

Wg{-)^lw,{-gY') + \w2{-g^') 

where Wi and W2 are the energies of the nearest-neighbour and next-nearest neighbour springs, respectively. 

But, when gi ^ g2, there is no g such that Wg{-) = Wi{-g^^). This is, of course, due to the interaction 
between the nearest-neighbour springs and the next-nearest-neighbour .springs. In dimensions larger than 
one, when shear resistance is present, such interaction is generic; indeed such interaction is the source of 
shear resistance. 

On the other hand, gi = g2, is analogous to isotropic growth. 

4.1.1. Comparison of the two decompositions for the lattice (Z'^, {ei, 62, ei±2}) 

To illustrate the difference between the energy-deformation decomposition (9) and the multiplicative 
decomposition (10) we consider the lattice (Z^, {ei, 62, ei±2}) with initial rest-lengths of the springs Li = 
L2 = 1 and Li+2 = -^1-2 = V^. Note that this is a residually unstressed configuration. 

In Examples 4.8 to 4.11 below we compute the fractional error 

that results when the energy of the recombined system is computed using the multiplicative decomposition 
for 

Ai A3 



Fe 



A. 



Ai, A2 e [0.8, 1.25], A3 e [-0.5, 0.5] } . (19) 



To do this, we should first associate a (single) G G M^^^ with the grown system. It is natural to identify 
G with the deformation gradient corresponding to the ground state of the grown system. This specifies G 
unto pre-multiplication by a rotation in SO{2). In keeping with the form of F in (19) we pick G to be of 
the form 

Gil G12 
G22, 



G 

and require Gn, G22 > 0. 



Example 4.8 (Recombination that shrinks only the non-diagonal springs by 10%). Let the recombination 
in the system be <7i = 52 = Ij 51+2 = 51-2 — 0.9. An elementary calculation shows that G = 0.94475/. 
(Note however that the growth is not isotropic.) The region where the fractional error (18) is greater than 
10% is shown in Figure 5a. 



Example 4.9 (Recombination that elongates only the diagonal springs by 10%). Let the recombination in 
the system be gi ~ 52 ~ 1, .91+2 = .91-2 = l-l- An elementary calculation shows that G = 1.04525/. (Note 
however that the growth is not isotropic.) The regions where the fractional error (18) is greater than 10% 
is shown in Figure 5b. 
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(a) (b) 




(c) (d) 

Figure 5: Regions where the fractional error in energy is greater than 10% for (a) Example 4.8, (b) Ex- 
ample 4.9, (c) Example 4.10 and (d) Example 4.11. Note that (c) is very similar to (d), showing the 
unimportance of residual stress; on the other hand these regions are very different from (a) and (c) showing 
the importance of shear caused by growth. This is in keeping with Remark 4.5. 
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Example 4.10 (Recombination of about 10% resulting in a residually- unstressed sheared lattice). Let the 

recombination in the system be gi = 92 = Ij .9i+2 = 0.9 and gi-2 = \/2 — 0.9^ ~ 1.09087. Note that 



(^51+2) +(\/2<?i-2) =2{gf+gl), 
so the recombined system is also residually unstressed. An elementary calculation shows that 



G = 



1 -0.19 
0.98178 



The region where the fractional error (18) is greater than 10% is shown in Figure 5c. 

Example 4.11 (Recombination resulting in a residually-stressed sheared lattice). Let the recombination in 
the system be gi = g2 = I7 31+2 — 0.9 and gi-2 = l-l- Note that 

[V2g,+2y + (\/2.gi_2)' ^ 2{gf+gi), 
so the recombined system is residually stressed. An elementary calculation shows that 

1.00243 -0.19757\ 
0.98277 J ■ 



The region where the fractional error (18) greater than 10% and 20% are shown in Figure 5d. 



4.2. Lattices on "Z"- . 

Next we turn to generic homogeneous lattices undergoing homogeneous growth. 

Consider the lattice (Z^,C). Let the initial rest-lengths of the springs he Ly, v £ C. Let the growth 
(due to recombination) in the springs he g^, v € C. The following result relates the continuum limits of the 
initial and grown systems: 

Theorem 4.12 (Homogenenously-recombined {Z^,C)). Let (Z^,C) be a lattice of order K and let its 
Cauchy-Born continuum energy density be Wi and Wg, before and after homogeneous recombination, re- 
spectively. Then there exist Wi, . . . , Wk ■ M.-^^-^ — > M (independent of ga, a £ C) and Gi, . . . , Gk G M-^^-^ 
(independent of La, a £ C) such that VF G R-^^-^, 

K 



W,{F) = J2Wk{F), (20a) 
fc=i 

K 

fc=i 

When K — 1, from Theorem 4.12 the energy-deformation decomposition reduces to a multiplicative 
decomposition of the deformation: 

Corollary 4.13 (Multiplicative decomposition of the deformation holds for a growing fluid). Let the as- 
sumptions of Theorem 4-12 hold for a lattice of order 1. Then there exists G E M.^^^ such thatVF G M^^^, 

WgiF) - W,{FG-'). 

Proof of Theorem 4.12. Let Ck, k — 1, . . . , K he & partition of C as specified in Definition 3.3. From (7), the 
Cauchy-Born continuum energy densities are of the form 

vec \ ^ / k=iveCk ^ ^ 

vec ^ '" ^ k=iveCk \ " / 
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Let Ck 3 Cfc be a basis of . We define Gk by 

G-^v^g;;^v, veCk, (22a) 

Gk^v = v, veCk\Ck. (22b) 

and Wk-. M^^-° -> M by 

Wk{F) -.^Y.^ C^) (23) 



veCk ^ 

Theorem 4.12 follows immediately from (21), (22) and (23) 



□ 



Remark 4.14. When D = 2 and C = {ei,e2,ei±2} we recover Theorem ^.i(l). It is clear that K ^2 and 
the possible partitions o/ {ei, 62, ei±2} o-^e 

{ei,e2},{ei±2}, 
{ei,ei+2},{e2,ei_2}, 
{ei, 61-2}, {62, 61+2} 

and the corresponding decompositions of the energy are given by (11), (12) and (13), respectively. 

5. Growth and homogenisation 

Next we explore how the energy-deformation decomposition interacts with homogenisation by considering 
an inhomogeneous lattice subjected to inhomogeneous recombination. The natural question is whether 
Theorem 4.12 extends to this situation, that is, are the continuum energies of the initial and recombined 
lattice related by an energy-deformation decomposition? 

Corollary 5.1 (of Theorem 4.12 below shows that indeed they are. We state it for inhomogeneous 
lattices with homogeneous representatives (Definition 3.5) because we do not wish to be distracted by the 
technicalities of homogenisation; a more general result is proved in Chencliiah and Sliipman (2013a). 

Corollary 5.1 (Recombination of inhomogenenous (Z^,C)). Let the initial and recombined (inhomoge- 
neous) lattices (Z^,C) lattices have homogeneous representatives whose Cauchy-Born continuum energy- 
densities are Wi and Wg respectively. Then there exist Wi,. . . , Wk ■ M.^^^ M and Gi, . . . , Gk G E-°^-° 
such that Vi^ G M^^^, 

K 

W',(F) = ^Tyfc(^^), 

K 

Wg{F) = J2Wk{FG-'). 

k=l 

Proof. By assumption (see Definition 3.5 and (7)) there exist Ci,Cg: C U — C — > (0, 00) such that Wi and 
Wg are the continuum energy densities of the homogeneous lattices (Z^,C, Ct) and {Z^,C, Cg), respectively. 
Thus we are in the situation of Theorem 4.12 and the result follows. □ 

In the rest of the section we present computational simulations that illustrate Corollary 5.1 for the 
lattice (Z^, {ei, 62, ei±2}) by demonstrating the existence of Wi, W2: M^^^ — > M and Gi,G2 G M^^^ such 
that VF e M2x2^ 

W,{F) = Wi{F) + W2{F), 
Wg{F) = Wi{FG;') + W2{FG~^) 
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(see (9)). Here Wi and Wg are the continuum energies of the initial and rehomogenised lattices. 

First, in Section 5.1, we consider lattices that are homogeneous prior to (but not after) recombination, 
and demonstrate that the energy-deformation decomposition is stable under homogenisation of growth. 
Then, in Section 5.2, we consider lattices that are inhomogeneous both before and after recombination and 
show that the energy-deformation decomposition is stable also under homogenisation of elasticity. 

5.1. Homogeneous {1,^1, {ei, €2, ei±2\) recombining inhomogeneously 

Deformation gradients. In the simulations below both the initial and the recombined system are subjected 
to afhne displacement at the boundary corresponding to a deformation gradient F (which is thus also the 
average deformation gradient in the system). This is picked from a set of the form 



Ai A3 
A2 



Ai,A2 € [^,A<i],A3 G hA„A,] (24) 



where Ad, A^ > 1. However, to permit two-dimensional plots the figures below are plotted for deformations 
of the form 

F\X)=(^^ J), Ae[^,A,], (25a) 

P'W-[o a)' ^^[x;^'^'^]' (25c) 
F'{X)^ (J 1)^ Ae[-A„A,]. (25d) 

(The superscripts stand for "horizontal", "vertical", "dilational" and "shear" respectively.) 

Existence ofWi,W2: E^^^ — >• E andGi,G2 G E^^^. Wi and 14^2 can be immediately picked from Lemma 4.3. 
For Gi,G2 S E^^^, we search for those that minimises the relative mean square error 



E 



' WgjF) ~ {WijFG^') + W2{FG^')) y 



(where the sum is approximated by a uniform sampling over T). We accept these to be the desired Gi, G2 
provided 

1. The error is less than a tolerance, and 

2. Gi and G2 are convergent, and the error is non-increasing, when the grid size increases. 

We also compare the resulting Gi, G2 with the best choices, for each A, for the deformations in (25). For 
example, G]^(A),G2(A) minimises 




WgiF^jX)) ~ (W-i(J^-(A)Gr^) + W2{F^{X)G^')) y 
WgiF^iX)) 



{26) 



and we compare it with Gi,G2. 

For the spring energy W (see (4b)) we pick 

W{-) = (• - 1)' (27) 

with q = 2,3,4, depending on the simulation. 
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-OJ -0.1 



D.1 0.2 



Figure 6: Results of Simulation 5.2. The true energy (black dots) and the homogenized energy (blue squares) 
are plotted as functions of A for deformations of the form (a) F'^(A) and (c) F''{X). The fractional error of 
the homogenized to true energy are plotted in (b) for the results in (a) and in (d) for the results in (c). 



Simulation 5.2 ((Zfg, {ei, 62, ei±2}) with alternating growth in diagonal springs). 

Initial state: Li — L2 — 1 and Li^2 = ^1-2 — V^. Note that the rest-state has zero energy, that is, is 
residual-stress free. 

Recombination: gi = §2 = 1; that is, the horizontal and vertical springs did not grow. 51+2 was chosen 
to alternate between 1.2 and \/2 — 1.2^ in a checkerboard fashion, and gi^2 was chosen to be y^2 — .9i_2- 
Note that each unit cell of the recombined lattice possesses a zero-energy state. 

Result: We choose W\ and W2 as in (11), with q — 2 in (27). Since gi — g2 — 1, we choose Gi to be 
of the form 71/. From (lib) G2 is of the form 

* V 71-2/ * 

Since the system possesses cubic symmetry we expect 71+2 = 7i-2; but we do not a priori impose this. 

Let us first impose a deformation F'^(\) for A G [1/1.25,1.25]: 7^ = 1.2010, 7^+3 = 7^-2 = 0.8300 
minimises the least-square relative error over all A in the chosen range. Figure 6(a) compares the true 
energy with the energy given by homogenisation for each X, and the fractional error of the homogenized 
energy to the true energy as a function of X is shown in Figure 6(c). 

Next we deform the system by F%X) for X G [-0.25,0.25]; 7^ = 0.8970, 7^+2 = 7i-2 = 0.8020 minimises 
the least-square relative error over all X in the chosen range. Figure 6(b ) compares the true energy with the 
energy given by homogenisation for each A, and the fractional error of the homogenized energy to the true 
energy as a function of X is shown in Figure 6(d). 

In Simulation 5.2, the tolerance associated with the choice of Gi, G2 is 4% error; in the simulations that 
follow the maximum error is about 10%. Table 1 gives computations of and 7i-|_2;7i_2 the mean 
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N 


if 


7l+2 — 7l-2 


mean squared error (26) for (Z^, {ei, 62, ei±2}) 


2^ 


1.2050 


0.8270 


0.1405 


24 


1.2010 


0.8300 


0.1676 


25 


1.1970 


0.8300 


0.1618 


26 


1.1980 


0.8300 


0.1600 



Table 1: Computations of jf and 7i+2:7i-2 fo^' simulations identical to Simulation 5.2 except for grid size. 
The mean squared error (26) is calculated over 60 equidistant choices of A between 0.8 and 1.25. 

squared error (26) for simulations similar to Simulation 5.2, but with grids (Z^, {ei, 62, ei±2}) of varying 
size. 

Remark 5.3. Simulation 5.2 is related to a situation that occurs in the growth of plant cells; see Baskin 
(2005). For a cluster of neighbouring cells with the same cellulose microfibril orientation, the growth of the 
plant cell is related to the microfibril orientation. However when there is a large variation in microfibril 
orientation in a cluster of cells, then microfibril orientation of a cell is less correlated with the growth of 
that cell. In this case, as in Simulation 5.2, the homogenised growth of the cell closer is related to the 
homogenised microfibril orientation. 

Simulation 5.4 ((Zjg, {ei, 62, ei-1-2}) with growth in all springs uniformly distributed about 1). 

Initial state: Li ^ L2 — 1 and Li-1-2 — L1-2 = V^- Note that the rest-state has zero energy, that is, is 
residual- stress free. 

Recombination: gi,g2,gi+2, and gi^2 were all chosen from uniform distributions on [0.8,1.2]. In 
contrast to Simulation 5.2, generically 

(^2.91+2)' + (V2gi_2)' ^ 2 {gl+gl) , 

so unit cells of the recombined lattice generically do not have a zero-energy state. Thus growth leads to 
residual stress. 

Result: We choose Wi and W2 as in (11), with q ^ 2 in (27). We choose Gi to be of the form 71/ and 
G2 to be of the form 72/. 

Let us first impose a deformation F'^{X) for A G [1/1.5, 1.5]; 7^ = 1.0240, 72 = 0.8600 minimises the 
least- square relative error over all A in the chosen range. Figure 7(a) compares the true energy with the 
energy given by homogenisation for each A, and the fractional error of the homogenized energy to the true 
energy as a function of X is shown in Figure 7(c). 

Next we deform the system by F'^{X) for A G [—0.5,0.5].' 7^ — 1.050 and 7I — 0.993 minimises the 
least-square relative error over all A in the chosen range. Figure 7(b) compares the true energy with the 
energy given by homogenisation for each A , and the fractional error of the homogenized energy to the true 
energy as a function of X is shown in Figure 7(d). 

Repeating this simulation, but choosing 51,52 from uniform distributions on [1 — Sg^v,^ + Sg^v] and 
5i+2j5i-2 from uniform distributions on [1 — Sgii, 1 + Sgd], we plot the resulting homogenization factors 7^ 
and 72 in Figure 8. 

Results using g = 3 (7^ = 1.0100, 7I = 0.9590/ q = 4 (^f = 1.0060, 7^ = 0.9640/ and deformations 
F'^ are shown in Fig. 9. 

Simulation 5.5 ((Z^g, {ei, 62, ei±2}) with residual stress; growth in all springs uniformly distributed about 
!)• 

Initial state: Li = L2 — I and ii+2 — ii-2 — 1.25. Note that the rest-state has non-zero energy, that 
is, it has residual stress. 

Recombination: gi, g2, gi+2, and g 1-2 were all chosen from uniform distributions on [0.8,1.2]. Since, 
generically, 

(1.25.91+2)' + (1.25gi_2)' ^ 2 (5^ + gl) , 
20 



(a) 



(c) 



1 1.1 1J 

A 



1 , 1.1 1.2 

A 




Figure 7: Results of Simulation 5.4 for q = 2. The true energy (black dots) and the homogenized energy 
(blue squares) are plotted as functions of A for deformations of the form (a) F'^{X) and (c) F''{1). The 
fractional error of the homogenized to true energy are plotted in (b) for the results in (a) and in (d) for the 
results in (c). 
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Figure 9: Results of Simulation 5.4 for (a,b) q = 3 and (c,d) q = 4. The true energy (black dots) and the 
homogenized energy (blue squares) are plotted as functions of A for deformations of the form F'^{X). The 
fractional error of the homogenized to true energy are plotted in (b) for the results in (a) and in (d) for the 
results in (c). 
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Figure 10: Results of Simulation 5.5. The true energy(black dots) and the homogenized energy (blue squares) 
are plotted as functions of A for deformations of the form (a) F'^{X) and (c) F''{1). The fractional error of 
the homogenized to true energy are plotted in (b) for the results in (a) and in (d) for the results in (c). 



generically, unit cells of the recomhined lattice do not have a zero- energy state. Thus the grown system is 
also residually stressed. 

Result: We choose Wi and W2 as in (11), with q — 2 in (27). We choose Gi to be of the form 71/ and 
G2 to he of the form 72/. 

Let us first impose a deformation F'^{\) for A G [1/1.5,1.5]; = 1.132, 72 = 0.957 minimises the 
least- square relative error over all A in the chosen range. Figure 10(a) compares the true energy with the 
energy given by homogenisation for each A, and the fractional error of the homogenized energy to the true 
energy as a function of X is shown in Fig. 10(c). 

Next we deform the system by F'^iX) for A G [—0.5,0.5]; 7f — 1.007 and 7I — 1.085 minimises the least- 
square relative error over all X in the chosen range. Fig. 10(b) compares the true energy with the energy 
given by homogenisation for each X , and the fractional error of the homogenized energy to the true energy 
as a function of X is shown in Fig. 10(d). 



5.2. Inhomogenenous (Z^, {ei, 62, ei±2}) recombining inhomogenenously 

In the simulations of Section 5.1 we were guided by Lemma 4.3 in our choice of the the energies Wi and 
W2 in the decomposition. For all three simulations we picked these as in (11), but (12) and (13) would have 
served just as well. 

When the initial lattice is inhomogeneous we determine Wi and W2 from the homogeneous representative 
of the initial lattice. Such a homogeneous lattice has a continuum energy of the form 

ie{l,2,H-2,l-2} 



La 



23 



and thus it remains to find the optimal homogenised rest-lengths La, a E {1,2,1 + 2,1 — 2}. We do this by 
picking these to minimise the relative mean square error 





2 



(where the outer sum is approximated by a uniform sampling over J-). We accept these to be the desired 
rest-lengths provided that 

1. The error is less than a tolerance, and 

2. the rest-lengths are convergent, and the error is non-increasing, when the grid size increases. 

We now allow the lattice to recombine (inhomogeneously) and proceed as in Section 5.1. 

Simulation 5.6 ((Z^g, {ei, 62, ei±2}) with nonuniform rest lengths; growth in all springs uniformly dis- 
tributed about 1). 

Initial state: Li and Li were chosen from uniform distributions on [0.8,1.2]. L1+2 and Li_2 were 
chosen from uniform distributions on [0.8-^/2, 1.2-\/2] • Note that the rest-state has non-zero energy, that is, 
it has residual stress. 

Recombination: 32, 51+2, and gi-2 were all chosen from uniform distributions on [0.8,1.2]. Since, 
generically, 

(1.25.91+2)' + (1.25gi_2)' ^ 2 (5^ + gl) , 

generically, unit cells of the recombined lattice do not have a zero-energy state. Thus the grown system is 
also residually stressed. 

Result: We choose W\ and W2 as in (11), with q = 2 in (27). 

Let us first impose deformations F'^{X) for A S [1/1.5, 1.5]; Homogenised rest lengths Li = L2 = 1.0820 
and L14.2 = Li^2 — 0.8790 minimise the least-square relative error over all A in the chosen range (with 4% 
maximum fractional error). Allowing the springs to grow, homogenised growths of the form Gi — 71/ and 
G2 = I2I with ~ 1.0560 and 72 = 1.0330 minimise the least-square relative error over all A in the chosen 
domain (with 5% maximum fractional error). 

Next we deform the system by F^{X) for A G [—0.5, 0.5]; Homogenised rest lengths Li ~ L2 ~ 1.0260 and 
Li_l_2 = -^^1-2 = 0.8530 minimise the least-square relative error over all A in the chosen range (with 2.5% 
maximum fractional error). Allowing the springs to grow, homogenised growths of the form Gi = 71/ and 
G2 = I2I with = 1.1930 and 72 = 0.8110 minimise the least-square relative error over all A in the chosen 
domain (with 6% maximum fractional error). 



6. Outlook: Recombining networks and continua 

While the precise analysis above is limited to lattices, the central insights regarding multiplicative de- 
compositions and energy-deformation decompositions hold for networks in general. (By a network we mean 
a system of springs connected through nodes, a lattice being a special case.) 

This follows from the observation that in a I?-dimensional network that is able to resist shears a generic 
node should be connected to at least D-\-l springs, see Figure 11. It follows, from the reasoning in Section 4, 
that the continuum energy of such a network cannot obey a multiplicative decomposition. 

On the other hand when each node of a network is connected to no more than 2D other nodes then, 
reasoning as in the proof of Theorem 4.12, we are led to hope that two energy-densities and two growth 
descriptors would describe the continuum energy-density of the network. Thus we hypothesise: 

Hypothesis 6.1 (D-dimensional shear-resisting networks with connectivity no more than 2D). The initial 
continuum energy- density Wi of a D-dimensional shear-resisting network with connectivity no more than 
2D is related to its recombined continuum energy- density Wg through 

W,{-)^Wi{-) + W2{-), (9a) 
Wg{-)^Wi{-G^^) + W2{-G:^^). (9b) 
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(a) 



(b) 



Figure 11: Shear-resisting networks in two and three dimensions where each node is minimaUy connected. 



for some Wi, W2 : M^^^ — > M, each of which is independent of growth and vanishes on a shear, and some 
Gi, G2 G M.^^^ which are independent of the deformation. 



From (9), 



It is convenient to set 



This gives, 

(Alternatively, from (9), 



Wg{F) = W,{FG^') + {Wi{FG^^) - Wi{FG^')) 
H := 

Tjf 

11 . — Lt2^i 5 

W'{F,-) := Wi{F-) ~ Wi{F). 
Wg{F) = W,{FH) + W'{FH,H'). 



(28) 
(29) 



Wg{F) = W,{FG^^) + {W2{FG^^) - W2{FG^^)) . 
and with the alternate definitions 

H := Gi^, 

Tjf 

11 . — LtiLt2 j 

W'{F,-) := W2{F-) - W2{F). 

we again obtain (29).) 

Observe that H and H' are measures of growth with the later measuring the extent to which growth is 
anisotropic. W' measures the extent to which the energy of the grown body differs from a motion of the 
initial energy. When the growth is isotropic then H' — I and W'{-, I) — thus recovering the multiplicative 
decomposition as a special case. 

We investigate this formulation further in Chcnchiah and Shipman (2013b). There we also investigate 
the question of which of the energy densities that have been proposed for soft tissues (see, for example, 
Bogen and McMahon (1979), Chaplain (1993), Delfino et al. (1997), Holzapfel et al. (2000), Sacks (2000), 
Horgan and Saccomandi (2002, 2003), Shergold et al. (2006)) are compatible with Hypothesis 6.1. One 
natural way to pursue this question is via the following definition: 



Definition 6.2 (Shear decomposability). W: 



fDxD 



w,. 



t,DxD 



is shear decomposable of order K if there exist 



— > M, i — 1, . . . , K , such that 
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1. Each Wi, i = 1, . . . , K , vanishes on a shear, and 

When K — 2 we shall omit "of order two" and say "shear decomposable". 

The significance of this definition is, of course, the following corollary of Hypothesis 6.1: 

Corollary 6.3. Let W : MP'^^ ^ M fee the energy- density of a recombining D-dimensional shear-resisting 
network with connectivity no more than 2D. Then W is shear decomposable. 

So the question can be rephrased: 

Question 6.4. Which of the energy densities that have been proposed for soft tissues are shear decomposable? 
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